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Abstract 

The construction of heavy quark effective field theory (HqEFT) is extended 
to arbitrary order in both expansion parameters a s and l/m q . Matching 
conditions are discussed for the general case, and it is verified that this 
approach correctly reproduces the infrared behaviour of full QCD. Choos- 
ing a renormalization scheme in the full theory fixes the renormalization 
scheme in the effective theory except for the scale of the heavy quark 
field. Explicit formulae are given for the effective Lagrangian, and one- 
loop matching renormalization constants are computed for the operators 
of order 1/m. Finally, the multiparticle sector of HqEFT is considered. 



1 Introduction 



The Heavy Quark Limit (HQL) [T[ has become a very useful tool for the de- 
scription of systems involving one heavy quark. The main progress consists in 
the exploitation of additional symmetries occuring in the limit of infinite mass 
of heavy quarks. These additional symmetries allow model independent predic- 
tions for systems with heavy quarks, e.g., they yield model independent relations 
between form factors of weak transition matrix elements and also model inde- 
pendent absolute and relative normalizations of form factors at certain kinematic 
points. 

The predictions obtained in the HQL receive corrections which are governed 
by two small parameters. The first type of corrections are the QCD short distance 
corrections which are obtained in a perturbation series in the parameter a s (m), 
the strong coupling constant at the scale of the mass m of the heavy quark. The 
second type are the recoil corrections governed by the small parameter Aqcd / fa. 

A convenient tool to deal with these corrections to the HQL is the so called 
Heavy Quark Effective Field Theory (HqEFT) which was originally formulated 
in [0, |3|, |J. The construction of any effective theory consists of two steps 0. The 
first step is to identify the degrees of freedom which are irrelevant at the scales 
under consideration and to remove them. In the language of functional integrals 
these degrees of freedom are integrated out in the functional integral. Although 
generally the action of the full theory is the integral of a local Lagrangian, this first 
step will leave us with a nonlocal action functional. However, this nonlocality 
is connected to the large mass M of the irrelevant degrees of freedom and an 
effective theory may be constructed, if the nonlocal action functional can be 
expanded into local terms where the expansion parameter is 1/M. 

This point of view has been elaborated in || for the case of HqEFT by identi- 
fying the massive degees of freedom and by explicitly integrating them out from 
the functional integral of QCD. In this way the tree-level contributions in all 
orders of 1/m have been obtained. In other words, this approach yields the cor- 
rect HqEFT at the scale of the heavy quark mass m, i.e., the tree-level matching 
of full QCD to the effective theory. However, computing the coefficients in the 
HqEFT Lagrangian to higher order in the perturbation expansion necessarily in- 
volves loop calculations using the Feynman rules of HqEFT. In the past these 
matching contributions usually have been obtained by evaluating corresponding 
diagrams in the full and effective theory separately, which is inconvenient because 
of spurious infrared (IR) divergences that cancel only in the final result. 

The second problem with this approach is connected to the multiparticle 
states of HqEFT. In the HQL the numbers of heavy quarks and heavy antiquarks 
are separately conserved and the derivation given in |J only deals with the one- 
particle sector of HqEFT. However, it has been noticed that in the two-particle 
sector some unusual features of HqEFT appear which are related to the fact that 
the naively calculated anomalous dimensions pick up an imaginary part which 
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leads to phases in the Wilson coefficient functions. It has been shown later p]] that 
one may in fact define real anomalous dimensions by redefining the multiparticle 
states of HqEFT in an appropriate way. 

The purpose of the present paper is to clarify these two points and to supple- 
ment the discussion given in ||. We shall extend the work in by stating the 
matching and renormalization conditions for the effective theory beyond tree level 
and thus provide a method to construct the effective theory to arbitrary order in 
the two expansion parameters. The ideas of |2| are extended in order to show that 
this theory indeed reproduces the results of full QCD in the given approximation. 
In Sec. 2 we fix our notation by showing at the level of Greens functions that the 
Lagrangian given in || correctly reproduces the QCD Greens functions at tree 
level. In Sec. 3 we include also loop effects and discuss the matching of HqEFT 
to full QCD. It is possible to summarize all corrections in an effective Lagrangian 
which takes a simple form, and thus to provide explicit formulae for the matching 
coefficients to all orders. As an application, in Sec. 4 we calculate the one-loop 
matching coefficients of the operators in the Lagrangian up to order 1/m. Finally 
we address the problems of multiparticle states in HqEFT and conclude. 



2 HqEFT at Tree Level 

We shortly review the derivation of the HqEFT at tree level as it has been given 
in |(J. The heavy quark fields occuring in the full QCD Lagrangian 

C = ^{iJp — m)ifj + fjip + ^yq (1) 

are rewritten using the projections on upper and lower components with respect 
to the reference frame given by the velocity vector v 

K = ^ (2) 

as 

P+tfj(x) = e- imv - x h v (x), (3) 

P~i/>(x) = e-"H v (x). (4) 

With the corresponding parameterization of the heavy quark sources 

p v (x) = P+e~77(x), (5) 
R v ( x ) = p- e imv - x T](x) (6) 

and after integration over the H v fields in the functional integral, the tree-level 
Lagrangian becomes a nonlocal expression 

4 0) = h v i{v D)h v + (h v iD^ + R v )- 1 ~{iD^h v + R v ) 

2m + i{v ■ D) — xe 

+ p v h v + h v p v , (7) 
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where the transverse derivative is given by 

D^ = ip-i,{vD). (8) 
The propagator of the heavy quark field can be read off as 

S+(k) = % ^ , (9) 
v - k + ie v ' 

and contains only forward propagation in time. 

The nonlocal effective Lagrangian may be expanded in orders of 1/m as 

4°) = h v i(v ■ D)h v + {KiJp + R v )—P^ £ i{iIph v + R v ) 

+ p v h v + h v p v , (10) 

where the backward propagation is lost if the series is truncated at any finite 
order in 1/m. In these expressions we have retained the sources R v of the lower 
component fields. They may be dropped if only Greens functions sandwiched 
between P+ projectors are to be calculated, which is true if the external states 
are eigenstates of the lowest order HqEFT Lagrangian. However, if the external 
state contains a P~ projection, differentiations with respect to R v are necessary. 

It may be verified that the Lagrangian ([IT]) correctly reproduces the tree 
level Greens functions of the full QCD Lagrangian to arbitrary order in the 1/m 
expansion. To see this, we rewrite the geometric series in (|T0|) as 



—i 00 / iv ■ D \ n 

E (~2^r) P v = S- + S-(ig4)S- + S-(ig4)S-(ig4)S- + .. . , (11) 

where the expression S v which becomes in momentum space 

S -f k \ = _ iP -_L f'UljA = ~ %P v (12) 
" 2m ~^ \ 2m J 2m + v • k 

may be interpreted as the propagator of the lower component field. To any finite 
order in the 1/m expansion this is a local expression and therefore it is part of 
the interaction terms in the effective Lagrangian. 

Using the Feynman rules as derived from (|TUD, with the upper component 
part of the Lagrangian in the form 

h v i{v ■ D)h v = h v P+i(v ■ d)P+h v + h v P v + (g4)P+h v , (13) 

where we have made explicit the projection operators, we may collect all terms 
that can occur in a Feynman diagram between two insertions of the gluon field 
gjf>.. If both vertices contain projectors, they sum up to 

— -f- (il/S v ift^J -\- ^i^S v ilfi) ^i^S v i^j -\- . . . 

ip+ l+v k/2m 
v vk + k 2 /2m' y ' 
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Similarly, if one insertion contains a P v projection, we have 

G^ = G^"^~ S v 

_ + ^/2m 

" v . k + k?/2m VJ [ ' 

whereas between two P~ projectors there can be an additional term from the 
geometric series ([11]) 

G~ = S-i^Gf + S; 

v-k/2m 

1 v vk + k 2 /2m' 1 } 

In the sum of the four possible contributions the projection operators disappear 

G~l + + G* + G v + + G v = i- — — ■ ^ o (17) 

and the propagator of the full theory is recovered. The same is true if one or 
both of the gluon field insertions are replaced by external sources. Note that in 
the latter case the term quadratic in R v becomes important. 

The 1/m expansion of a tree-level Greens function is unique. We have seen 
that the full propagator, and thus any Greens function, is correctly reproduced 
by the effective theory. Truncating the series in the effective Lagrangian at a 
given order l/m k , we get an approximation to the full Greens function which 
differs by terms of order l/m h+1 : 



4?(A,m) - G£ } (A, 1/m) = O (g h+l - 2 X\X/m) 



k+l 



where 5 is the mass dimension of the Greens function, and the momenta of the 
h external heavy and I external massless lines are parameterized by 

Pi = mv + Xxi, 

Pi = Xyi, (19) 

so that the parameter A measures the deviation of the external momenta from 
mass shell. 

Although these facts might seem trivial for the tree-level case, we note that 
some parts of the full theory propagator shrink to a point and become part of 
the interaction vertices in the effective theory. Thus, the one-particle irreducible 
(1PI) Greens functions of the full theory do not correspond to the 1PI functions 
in the effective theory. This fact has to be accounted for in the matching of higher 
dimensional operators at order a s or higher. 
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3 Matching and Renormalization Conditions 

If we want a statement like ([Tj|) to hold at higher orders in the loop expansion, 
we have to include in the effective Lagrangian all possible local operators up 
to a given mass dimension and heavy-quark number that cannot be excluded 
by symmetry arguments. Furthermore, the operators that are already present 
in the tree-level Lagrangian will be modified by multiplicative renormalizations 
which in general are ultraviolet (UV) divergent. However, no IR infinities can 
occur if the calculations are organized appropriately, so one can dispense of all 
IR regulators in the matching calculations. 

In comparing the one loop Greens functions of the full and effective theories 
we first consider the Greens functions which contain no external heavy quark. 
Diagrams without any internal heavy quark line are identical in both theories. 
On the other hand, diagrams with heavy quark loops have no counterpart in 
the effective theory. This is true for any diagram even if two heavy quarks with 
different velocities are present: The heavy quark propagator in coordinate space 



describes a particle moving along a classical straight world line in the forward 
light cone. If two heavy propagators are connected at one point (the origin) 
in coordinate space, the loop cannot be closed at any other point, so the loop 
integral vanishes. To be precise, because at the origin two distributions are 
multiplied, it is equal to some undefined constant which can be got rid of by 
a simple renormalization. In contrast to light particles, a heavy particle loop 
cannot introduce any nonlocal interaction. 

Thus, one can calculate the coefficients of operators without heavy quarks 
such as 



by computing the corresponding loop diagrams in the full theory. This is just 
the result we would have obtained if we had integrated out the heavy quark 
completely. At one loop the result is given by the determinant of the heavy 
fermions in the functional integral M. 

For Greens functions with external heavy quarks we need additional countert- 
erms. These are to be chosen in such a way that the condition analogous to (|TS|) 



for n loop Greens functions with h external heavy-quark and / external gluon 
linesQ 

A<?(A,ro) = Gfc\\,m) - Gff (A, 1/m) = O (^V;A 5 (A/m) fc+1 ln"(A/m)) , 



1 Throughout this paper we ignore the ghost fields, and do not write out the pure gluon 
and gauge fixing terms in the Lagrangian. The ghost fields can be treated in the matching 
in the same way as the gluon fields. To order 1/m the question of BRS invariance has been 
addressed in 0. A complete discussion to all orders is beyond the scope of this paper. 




(20) 




(21) 



(22) 
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can be satisfied for h < ho, if we know that it is already true in lower order of 
the loop expansion. For the lowest order (l/m°) Greens functions with h = 2 
this has been shown in 0. Extending the ideas developed there we shall show 
in the following paragraphs how the matching can be accomplished to arbitrary 
order in the 1/m expansion. 

For 5 + k > 0, ([22|) is equivalent to the requirement that the first S + k 
derivatives of A^(A,m) with respect to external momenta vanish for A — > 0. 
If the quantity A^(A,m), which is defined in the same way as A^(A,m) but 
includes counterterms only up to n — 1 loop order, admits an expansion of the 
form 

\S+k 

AffiX, m) = C m s + dXm 5 - 1 + . . . + C s+k ' 



(2m) k 

+ O (g h+l - 2 a n s \\\/m) k+l ln n (A/m)) , (23) 

this can be accomplished by introducing local counterterms of order a™ and up 
to order 1 / m k into the effective Lagrangian. Since u = 5 + k is the maximal UV 
degree of divergence of the effective theory diagrams, the coefficients Co to Cs+k 
which translate into coefficients of counterterms are in general UV divergent. To 



make (|22"D hold, we have to include counterterms exactly for those diagrams with 
uj > 0, so the set of matching conditions ( |2"2"D uniquely defines a renormalization 
scheme in the effective theory. 

If ( f23"|) did not hold, we would encounter IR divergences in evaluating Aff or 
its derivatives for vanishing external momenta. These could come from regions 
where one or more loop momenta become small. We consider a particular diagram 
of the full theory together with its counterparts in the effective theory up to 
order l/m k , including all counterterm insertions up to n — 1 loop order. Let 
us investigate the behaviour of the integrand in the region where the momenta 
of some subset of s massless propagators become small of order A. Since the s 
light propagators are the same in both theories, we can factor them out in the 
difference. The remaining expression corresponds to a Greens function — not 
necessarily connected — in lower order of the loop expansion. It is represented 
by the same set of diagrams, where the s light lines have been cut, their momenta 
being regarded as external. If the remaining n' loop momenta are integrated over, 
according to the induction hypothesis this Greens function satisfies 

A&;>(A,m) = O (g h+l - 2 a :'X s '(\/m) k+1 \n n \X/m)) , (24) 

where 8' = 5 + 2s — 4(n — n'). The total IR degree of divergence in this integration 
region is therefore 

p = -1(5' + k + 1) - 2s + 4(n - n')] = -5 - k - 1, (25) 
and thus no IR singularities show up in the first 5 + k derivatives of Aff . 
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Any diagram with 5 + k < automatically satisfies (E^) without additional 
counterterms, since by the same reasoning the matching conditions for its sub- 
diagrams exclude unexpected positive powers of m coming from UV divergent 
subintegrations. Simple power counting then directly gives (^2]) for these dia- 
grams. 

We did not consider additional IR divergences from integration regions where 
some momenta of heavy quark lines become small. In fact there are none. If we 
cut a heavy quark line, the remaining subdiagram has a larger number of external 
heavy quark lines and thus it is part of a different sector of the theory. Violations 
of the IR power counting could occur if this subdiagram is UV divergent, which 
can happen first at order 1/m 2 . We are tempted to renormalize it by imposing 
matching conditions also for Greens functions with more than h external heavy 
lines, but as explained in Sec. 5, their counterterms can contain divergent phases 
which should not appear in the Lagrangian. Fortunately, these terms give van- 
ishing contributions to the final result, because in evaluating them, we have to 
close the heavy quark line again and get a heavy quark loop which is zero by 
definition. We conclude that the multiparticle sectors are completely irrelevant 
for the renormalization and matching of the one-particle sector, and all possible 
IR singularities have been accounted for in the previous paragraphs. 

In practice, the counterterms needed in the effective Lagrangian are given 
by the 1/m expansion up to a certain order l/m k of the 1PI diagrams in the 
full theory, with all corresponding diagrams of the effective theory up to the 
same order subtracted. Since the effective theory does not contain any mass 
scale, the subtraction terms are defined uniquely by the requirement that the 
first 6 + k derivatives of the subtracted diagram are IR finite, where S is the mass 
dimension of the diagram, so that the effective theory is not needed explicitely 
in the matching calculation. 

In evaluating the counterterms introduced by a particular 1PI diagram into 
the effective theory, we first subtract all counterterm diagrams necessary to make 
this diagram UV finite in the full theory. These may be taken in the MS scheme 
as usual, with one exception: In order to have a renormalization group invariant 
mass as expansion parameter, we must also include a counterterm for the finite 
part of the quark mass renormalization. In this way the self-energy diagram 
contains no constant term, and the 1/m expansion is done in terms of the pole 
mass which has to be redefined in each order of the perturbation expansion. 
Furthermore, the counterterms in the effective Lagrangian all have nonpositive 
mass dimensionP]. 

We may put the results of the preceding paragraphs together in a formal 
language, and give an explicit expression for the effective theory counterterm of 
a particular diagram T with mass dimension 5(T). For a one-loop diagram, it is 

2 The possibility of expanding in terms of a different mass, which introduces a residual mass 
term into the effective Lagrangian, has been considered in 



8 



given by the integral of the expression 

Cr = Ts(r)+klr = (1 — J P )%{T)+kIvi (26) 

where Jr is the integrand of the full theory diagram. The operator %+k gives 
the Taylor expansion of the integrand in terms of external momenta up to the 
required order 5 + k. The operator J p gives the first p + 1 terms of a Laurent 
expansion in terms of internal momenta, where p is the IR degree of divergence of 
the term on which J p acts. This takes into account the diagrams of the effective 
theory. All expansions are done about zero momentum k = 0, where the heavy 
quark momenta are parameterized by p = mv + k. 

Because of the additional subtractions, the subtracted integrand Cr is no 
longer UV finite, but has an UV degree of divergence u = 5 + k. Since it is a 
polynomial in external momenta of order 5 + k which is designed to cancel the 
leading terms in the difference of the full and effective theory diagrams, it also 
exactly cancels the UV divergences in this difference. 

Evaluating the counterterms for n-loop graphs, we have to take into account 
all lower order counterterms in the effective theory. We observe that a recursive 
procedure emerges closely analogous to the well-known renormalization of UV 
divergences. Thus we may give the effective theory counterterm of any n-loop 
1PI diagram in closed form as an extension of Zimmermann's forest formula |]10| 



Cr - %{v)+k E II (~^(7)+fc) ^r, (27) 

F: FgFjeF 

where the modified Taylor operator T is defined as in (^), and the sum runs 
over all forests F of subdiagrams 7 not containing T itself. The "renormalized" 
integrand 

#r = E II (-^( 7)+fe ) It (28) 

all F 7S-F 

is the full theory diagram with all effective theory diagrams subtracted, i.e., a 
term contributing to the difference A$ which obeys the matching condition 
(|22|). Of course, we can apply these formulae also to the h = case, where the 
IR subtractions first appear in two-loop order. 

This treatment of the matching conditions is not just formal. In fact, it is 
convenient because all integrations are manifestly IR finite, and the integrals can 
be evaluated for vanishing external momenta. In dimensional regularization, the 
l/e n poles are due only to UV singularities. 

In the effective theory we have the freedom left to rescale the heavy quark 
field by a finite amount z. The matching condition ( p2| ) is then modified to 



Gff (A,m) - z h ^(X,l/m) = O (g h+l - 2 ^\\\/m) k+l W\\/m)) . (29) 



If the renormalization scheme of the full theory is fixed, choosing a value for 
z and requiring that (53) holds uniquely defines a renormalization scheme for 
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the effective theory. Using the MS scheme in the lowest order effective theory 
amounts to specifying a wave function renormalization constant z that is given 
by the finite part of the difference of the heavy quark self energy diagrams in the 
full and effective theory. In higher order of the 1/m expansion we have to insert 
the counterterms as given by (^7|) explicitely, with the self-energy contribution 
of the external lines factored out. 

We may collect all integrands Cr (|27| ) that contribute to the counterterm of 
some 1PI Greens function with h heavy-quark and I gluon legs. After integrating 
over internal momenta, this becomes a function Chi which is a polynomial of order 
5 + k in external momenta. Fourier transforming back into coordinate space and 
inserting gluon fields (we do not write out colour indices), we define the quantity 

C/i = E^c«"" w (^)An--^« (3°) 

i L 

which looks like a generating functional of 1PI vertices, but does not contain the 
heavy quark fields yet. In contrast to the effective action of the full theory, it 
can be expanded as a series of powers of d/m. By definition, we also include the 
tree-level contribution into C2 which is just the ordinary vertex gjji. 

In the following paragraphs, we consider only the case h < 2, the one-particle 
case. With the caveat discussed at the end of Sec. 5, the generalization to mul- 
tiparticle Greens functions is straightforward. 

The effective vertices which are summarized in C2 can be linked together by 
propagator projections in the same way as the vertices gA 1 at tree-level. Repeat- 
ing the reasoning of Sec. 2, we find that we can incorporate all loop effects in the 
effective Lagrangian by defining a generalized covariant derivative 

if = iQ + C 2 = tip + 0(a s ) (31) 

with projections 

iV+ = P v + ilpP v + = (ivD)P v + + 0(a s ), 

iV^ = P v + ifp- + P v - i fP+ = i D^ + 0(a s ), 

iV- = p-ifp- = -(wD)P v - + 0(a s ). (32) 

The effective Lagrangian which generalizes (0) then reads 

C v = h v iV+h v + (h v iV^ + R V )- ——(iT^hv + Rj+Co 

2m — %D V 

+ p v h v + h v p v , (33) 

where Co summarizes the counterterms without heavy quarks. This effective 
Lagrangian is valid to arbitrary order in the 1/m and loop expansions. 
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It is straightforward to apply this procedure to operator insertions. For in- 
stance, a heavy-light current K = qTijj, where q is some light quark field, is 
matched onto 

' (W^h v + Rv) + h v 



K v = qJC 



2m — iDz 



(34) 



where /C stands for the sum of all 1PI diagrams of the full theory that involve one 
insertion of K, with an arbitrary number of gluons, where the IR subtractions as 
in (^) have been carried out. 

The form of C v and K v is restricted if BRS invariance holds in the matching. 
One would expect the generalized covariant derivative Tp to depend on the gluon 
field only through the ordinary covariant derivative iJp = + gjfi, at least in the 
background field gauge JT^] . This has been shown for the lowest order (l/m°) 
theory in ||. 

Furthermore, there are restrictions following from reparameterization invari- 
0. 
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They emerge from the fact that the Greens functions of the full theory 
do not depend on the heavy-quark velocity v and the residual momentum k sep- 
arately, but only in the combination p = mv + k. In the effective theory, the 
Lagrangian therefore depends on the velocity and the covariant derivative only 
in the combination V = v + iD/m. This introduces relations among operator 
coefficients in different orders of the 1/m expansion. 



4 Coefficients to order 1/m and a 



As an example, in this section we shall discuss the one-loop matching up to 
order 1/m. By combining denominators and algebraically reducing vector and 
tensor integrals, all quantities that occur in the one-loop matching calculations 
to arbitrary order in the 1/m expansion can be expressed in terms of IR-finite 
integrals in the form 

Ha6c(m) = / dk {k^Y{vmvk + k^/2mY' (35) 

of which only a few have to be evaluated because of the relations 

d c 
—E abc (m) = — — S a _i 6c+ i(m) (36) 
am 2m 2 

and 

Safe c (m) = 2m (H a+ i j6>c _i(ra) - 5 a+ i i6 _ liC (m)) . (37) 

To order 1/m there are three independent operators involving two heavy 
quarks possible which are usually written in the form 

J-h v {iD) 2 h v , ^-h v {iv ■ D) 2 h v , -£-h v a^G^h v . (38) 
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The second operator vanishes when sandwiched between physical states due to 
the heavy-quark equations of motion. 

Because of BRS invariance one only has to calculate the two one-loop vertex 
diagrams of the full theory, with the appropriate IR subtractions, to obtain the 
matching coefficients up to order a s and 1/m. Applying the formulae ( |26|j30|j33| ) , 
and expanding the result up to order 1/m, we obtain the one-loop effective 
Lagrangian at the matching scale \i — m in the background field gauge with 
£ = 1 (Feynman gauge) 



Hi) 



Z^p 1 Zhh v 



iv ■ D 



z m 2 _ z & ■ Df 



2m 



2m 



Urn M 



h, 



(39) 



where the matching renormalization constants are given by [1/e = 1/ e+ln(47r) /2- 
7_e/2, where the space-time dimension is 4 — e] 



Z h 
Zx 

z 2 



1 + 
1, 

1 



a* 



71 



O'.s 



7T 



+ 1 c. 



a, 

1 + — 

7T 



-)c / 



(40) 



and the wave-function renormalization of the full theory is in the MS scheme 



(41) 



The values of Z^, Z±, and Z 2 are also obtained by computing the one-loop self- 
energy to the required order. In addition, this gives the mass renormalization 



so that the 1 /m-expansion is done in terms of the pole mass 

m = m — 5m. 



(42) 



(43) 



The fact that Z\ — 1 is a consequence of reparameterization invariance [ |l3|j . 
BRS invariance ensures that no gauge-variant operators appear, and that only 
the unphysical constants Z^, Z h , and Z 2 depend on the gauge parameter. Our 
results for the one-loop matching renormalization constants agree with those 
given in WM. The calculational method of ||14|| , where dimensional regularization 



has been used for both UV and IR singularities, is thus justified by employing a 
scheme that is manifestly free of IR divergences. 

In order to obtain the effective Lagrangian valid at scales \i < m, one may sum 
up the leading mass logarithms by evaluating the UV-divergent parts of the loop 
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diagrams in the effective theory. The leading log calculation for the Lagrangian 
(|39D has been performed in |15|]. In fact, the one-loop anomalous dimensions are 
determined by the 1/e-poles of the renormalization constants (^0|). The finite 
terms in (fiTf) become relevant when subleading logarithms are summed. 



5 Hadronic States 

In a typical application of HqEFT one considers the matrix element of some 
operator O, describing e.g. the weak decay of a heavy quark, between hadronic 
states | in) and |out) involving heavy quarks. As in ordinary perturbation theory, 
the full QCD eigenstates are evolved from the corresponding lowest order HqEFT 
eigenstates |in(out)H) by adiabatically switching on the 1/m terms in the effective 
Lagrangian: 

(out|£>|in) = (out H| T [exp(-iS' v ) 6] |in H), (44) 

where S' v = J dx C' v {x) is given by the effective Lagrangian (|33|) excluding terms 
of order l/m°, and O is the result of a matching calculation such as fl54"D- 

Usually, fl44[) cannot be evaluated completely in terms of Feynman diagrams 
because perturbation theory breaks down at low scales. One rather uses per- 
turbation theory to sum up the leading (and subleading) logarithms \og(m/p), 
where p is some low hadronic scale, and expresses the result in terms of some 
operator O(p) which has the same matrix element between scaled-down states 
|in(out)H, p). These are defined in such a way that the matrix element is p 
independent: 

(out H| T [exp(-iS' v ) 6] |in H) = (out H, p\0(p)\in H, p). (45) 

Obviously, O(p) contains nonlocal operators (time-ordered products) in its 1/m 
expansion. In the solution of the renormalization group equations they mix with 
local operators. 

If we want to include into this treatment also states involving more than one 
heavy quark at a time, we have to extend the heavy-quark Lagrangian accord- 
ingly. We can apply HqEFT to asymptotic states where the heavy quarks end 
up in different hadrons, so each one can be assigned a separate velocity. The 
Lagrangian then contains a sum over all heavy quark and antiquark flavours and 
velocities that occur in the problem 

£ = E \h f JKH + (hlir$ + R f v ) 2m X _ w _ avM + Rf) 



f,v 



+ p(,hl + hlpl^j + E Co + (counterterms with h > 2), 



(46) 



where the projections of the generalized covariant derivative Tp (|31|) have been 
defined in (^). The counterterms without heavy quarks Co are inserted only 
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once for each quark flavour. The bilinear matching corrections that are incorpo- 
rated in Tp are needed once for each flavor and velocity. In addition, one expects 
counterterms involving more heavy quarks at higher order in the 1/m expansion 
if there is more than one heavy quark present in the initial or final states. These 
counterterms can link different velocity and flavour sectors. They are needed 
only if they involve no more heavy quarks than the process under considera- 
tion: Any diagram which contains an operator with more heavy fields than there 
are available as external lines necessarily contains a heavy quark loop and thus 
vanishes. 

If more than one heavy quark is present at the same time, some Feynman 
amplitudes develop UV divergent phases analogous to the well-known Coulomb 
phases of QED |ff6] . They can be absorbed into the definition of the multiparticle 
states itself [[?]]. In particular, the phase of the two-particle final state is related 
to the static interquark potential |7[]. If two velocities of different heavy quarks 
become equal, these phases become infinite, so that the whole expression is ill- 
defined. However, one could go back into coordinate space and consider the case 
of two separated heavy quarks which act as static coulour sources. This notion 
was taken in the first perturbative evaluation of the interquark potential |17|. 



6 Conclusions 

In the present paper we have extended the construction of HqEFT to arbitrary 
order in both expansion parameters 1/m and a s . Although our arguments as 
given in Sec. 2 are merely heuristic, we believe that they suffice to show that 
the effective theory reproduces the behaviour of full QCD near mass shell to any 
required accuracy. The formulae fl2~7f) and (^) provide a method to obtain the 
terms in the effective Lagrangian directly from the full theory. The fact that no IR 
regulators need to be introduced simplifies practical calculations. As an example, 
we have rederived the coefficients of the terms in the effective Lagrangian up to 
order 1/m and a s . 

For the construction of the complete effective Lagrangian ( |33| ) one needs in 
principle the knowledge of all diagrams of the full theory. One might argue that 
for this reason the effective theory is without physical content. However, its 
main power lies in the possibility of summing up large logarithms of the heavy 
quark mass, which are difficult to extract from the full theory expressions, and 
thus reorganizing the perturbation series. The method for calculating anomalous 
dimensions and solving the renormalization group equations within the effective 
theory is well known. The Lagrangian (|33| ) provides in closed form the initial 
conditions valid at the matching scale fi = m. It may also serve as a starting 
point for theoretical considerations, e.g., the proof of BRS invariance to all orders. 

In the last section we have considered the extension of HqEFT to sectors with 
more than one heavy quark. It should have become clear that these sectors do not 
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introduce difficulties if the heavy quarks remain well separated in the asymptotic 
states, because the UV-divergent phases that appear in these sectors can then 
be absorbed in the definition of the multiparticle states. 

All arguments in the present paper have been in the context of perturbation 
theory. We did not adress the question whether a nonperturbative approach to 
the 1/m expansion such as lattice gauge theory introduces any new problems. 
However, the perturbative matching of the continuum theory to the lattice reg- 
ulated version should proceed mainly along the same lines as we have discussed. 
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